from scipy.stats import poisson

# 定义泊松分布的参数 λ（平均发生率）
lambda_ = 15  # 平均每天15位顾客

# 定义我们感兴趣的事件发生次数 k
k = 20  # 我们想要计算某一天恰好有20位顾客的概率

# 使用poisson.pmf方法来计算概率质量函数 P(X = k)
probability = poisson.pmf(k, lambda_)

# 打印结果
print(f"在平均每天{lambda_}位顾客的情况下，某一天恰好有{k}位顾客的概率是 {probability:.6f}")